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SPECIFICATION 



Electronic Version 1.2.8 
Stylesheet Version 1.0 

METHOD AND APPARATUS FOR 
RECONSTRUCTION OF IMAGES 
IN PARALLEL MRI SYSTEMS 

Background of Invention 

[0001] The invention relates generally to Magnetic Resonance Imaging (MRI) systems and 
more specifically to a method and apparatus for reconstructing images in parallel MRI 
systems. 

[0002] Generally, MRI is a well-known imaging technique. A conventional MRI device 
establishes a homogenous magnetic field, for example, along an axis of a person's 
body that is to undergo MRI. This homogeneous magnetic field conditions the interior 
of the person's body for imaging by aligning the nuclear spins (in atoms and 
molecules forming the body tissue) along the axis of the magnetic field. If the 
orientation of the nuclear spins is perturbed out of alignment with the magnetic field, 
the nuclei attempt to realign their spins with the axis of the magnetic field. 
Perturbation of the orientation of nuclear spins may be caused by application of radio 
frequency (RF) pulses tuned to the Larmor frequency. During the realignment process, 
the nuclei precess about the axis of the magnetic field and emit electromagnetic 
signals that may be detected by one or more RF detector coils placed on or about the 
person. 

[0003] f requenC y 0 f the magnetic resonance (MR) signal emitted by a given 

precessing nucleus depends on the strength of the magnetic field at the nucleus 1 
location. As is well known in the art, it is possible to distinguish signals originating 
from different locations within the person's body by applying a gradient to the 
magnetic field across the person's body. For the sake of convenience, direction of this 
field gradient may be referred to as the left-to-right direction. Signals of a particular 
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frequency acquired during application of the field gradient may be assumed to 
originate at a given position within the field gradient, and hence at a given left-to- 
right position within the person's body. The application of such a field gradient is also 
referred to as frequency encoding. 

[0004] However, the application of a field gradient does not allow for two-dimensional 

spatial encoding, since all nuclei at a given left-to-right position experience the same 
field strength, and hence emit signals of the same frequency. Accordingly, the 
application of a frequency-encoding gradient, by itself, does not make it possible to 
discern signals originating from the top versus signals originating from the bottom of 
the person at a given left-to-right position. Spatial encoding has been found to be 
possible in this second direction by application of gradients of varied strength in a 
perpendicular direction prior to acquisition of the signal, to thereby twist the phase of 
the nuclear spins by varied amounts. The application of such additional gradients is 
referred to as phase encoding. 

[0005] Frequency-encoded data sensed by the RF detector coils after a phase encoding 
step is stored as a line of data in a data matrix known as the k-space matrix. Multiple 
phase encoding steps are performed in order to fill the multiple lines of the k-space 
matrix. An image may be generated from this matrix by performing a two- 
dimensional Fourier transformation of the matrix to convert this frequency 
information to spatial information representing the distribution of nuclear spins or 
density of nuclei of the image material. 

[0006] Alternatively, spatial encoding can be performed in three dimensions by applying 
phase-encoding gradients in two orthogonal directions, followed by a frequency- 
encoding gradient in the third orthogonal direction, with signals acquired during the 
frequency-encoding gradient, in order to generate a three-dimensional matrix of k- 
space data. Three-dimensional Fourier transformation of the matrix will then convert 
this frequency information to spatial information representing the distribution of 
nuclear spins or density of nuclei within a volume of the image material. 

[0007] Imaging time is largely a function of desired signal-to-noise ratio (SNR) and the 
speed with which the MRI device can fill the k-space matrix. In conventional MRI, the 
k-space matrix is filled one line at a time. Although many improvements have been 
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made in this general area, the speed with which the k-space matrix may be filled is 
limited. To overcome these inherent limits, several techniques have been developed to 
effectively simultaneously acquire multiple lines of data for each application of a 
magnetic field gradient. These techniques, which may collectively be characterized as 
"parallel imaging techniques", use spatial information from arrays of RF detector coils 
to substitute for the encoding which would otherwise have to be obtained in a 
sequential fashion using field gradients alone. The use of multiple detectors has been 
shown to multiply imaging speed, without increasing gradient switching rates or RF 
power deposition. Parallel imaging techniques fall into one of two categories. They 
can fill in the omitted k-space lines prior to Fourier transformation, by constructing a 
weighted combination of neighboring lines acquired by the different RF detector coils. 
Or, they can first Fourier transform the limited k-space data to produce an aliased 
image from each coil, and then unfold the aliased signals by a linear transformation of 
the superimposed pixel values. In either case, the reconstructed images tend to suffer 
from incomplete removal of aliasing artifacts, especially for large speedup factors. In 
images corrupted by aliasing, the edges of the image are wrapped into the center of 
the image. 

[0008] Another problem in the reconstruction of the image is degraded signal-to-noise 
ratio (SNR) in some regions of the image. One reason for both of the above effects is 
the overlapping of the sensitivity profiles of the coils in the array. Insufficiently 
distinct sensitivity profiles can lead to both residual aliasing and degraded regional 
SNR. One way to overcome the above problems is to keep the sensitivity profile of 
each coil orthogonal to all of the others. 

[0009] It is therefore desirable to reconstruct an artifact-free image, by reducing the 
effects of aliasing and degraded SNR. 

Summary of Invention 

[0010] Briefly, in accordance with one embodiment of the invention, a reconstruction 
method for use in a parallel MRI system is provided. The MRI system comprises a 
plurality of MR detector coils arranged in an array and each coil has a corresponding 
spatial sensitivity profile. The method comprises detecting a plurality of MR signals 
from the plurality of MR coils and processing the detected MR signals with at least one 
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filter bank to reconstruct at least one image. 

[001 1] In an alternate embodiment, a parallel Magnetic Resonance Imaging (MRI) system 
is provided for reconstruction of images. The MRI system comprises an array of 
magnetic resonance (MR) detector coils arranged in an array for detecting a plurality 
of MR signals, each of the coils having a corresponding spatial sensitivity profile and a 
processing means for processing the plurality of MR signals with at least one filter 
bank to reconstruct at least one image. 

Brief Description of Drawings 

[0012] These and other features, aspects, and advantages of the present invention will 
become better understood when the following detailed description is read with 
reference to the accompanying drawings in which like characters represent like parts 
throughout the drawings, wherein: 

[001 3] FIG. 1 is a schematic block diagram of an exemplary MR imaging system suitable 
for use with the present invention embodiments ; 

[0014] FIG. 2 is a block diagram illustrating a method for reconstructing images by a MRI 
system ; 

[001 5] FIG. 3 is an exemplary embodiment of a lapped transform filter bank used for 

reconstruction of images when the sensitivity profiles of the array elements overlap; 

[001 6] FIG. 4 is a schematic block diagram of an exemplary embodiment for a strip array; 

[001 7] FIG. 5 is a block diagram of a mathematical equivalent of the embodiment shown 
in FIG. 4; and, 

[001 8] FIG. 6 is an exemplary block diagram of a reconstruction for the embodiments 
shown in FIGURES 4 and 5. 

Detailed Description 

[0019] 

FIG. 1 is an example embodiment of a magnetic resonance (MR) imaging system 
1 00. MRI system 1 00 includes a computer 1 1 0, which controls gradient coil power 
amplifiers 1 30 through a pulse control module 1 20. The pulse control module 1 20 
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and the gradient amplifiers 1 30 together produce the proper gradient waveforms Gx, 
Cy, and Cz, for either a spin echo, a gradient recalled echo pulse sequence, a fast spin 
echo, or other type of pulse sequences. The gradient waveforms are connected to 
gradient coils 16, which are positioned around the bore of an MR magnet assembly 
130 so that gradients Gx, Gy, and Gz are impressed along their respective axes on the 
polarizing magnetic field B from magnet assembly 140. 

[0020] The pulse control module 1 20 also controls a radio frequency synthesizer 1 78 
that is part of RF transceiver system 1 70, portions of which are enclosed by dashed 
line block. The pulse control module 120 also controls an RF modulator 179, which 
modulates the output of the radio frequency synthesizer 1 78. The resultant RF 
signals, amplified by RF power amplifier 1 80 and applied to a RF coil assembly 1 47 
through transmit/receive switch 160, are used to excite the nuclear spins of the 
imaged object (not shown). 

[0021] In embodiments of the present invention, RFcoil assembly 147 is a RF detector 

array consisting of multiple detector elements, such as loops or conductive strips and 
may also be configured as an array of multiple coils or conductive strips. Other 
conductive materials and structures, for example copper rods, pipes, wires or other 
similar line structures, may also be used as detector elements. In a further 
embodiment, RF coil assembly 147 is an MRI phased array. In an embodiment, the 
detector elements may be non-overlapping or overlapping. 

[0022] The MR signals from the excited nuclei of the imaged object are picked up by RF 
coil assembly 147 and presented to preamplifier 171 through transmit/receive switch 
160, to be amplified and then processed by a quadrature phase detector 173. The 
detected signals are digitized by a high speed A/D converter 1 75 and applied to 
computer 1 10 for processing to produce MR images of the object. Computer 1 10 also 
controls shimming coil power supplies 141 to power shimming coil assembly 143. 
The manner in which computer 1 10 processes the detected signal to reconstruct an 
image is described in further detail below. 

[0023] Fig. 2 is a block diagram illustrating a method for reconstructing images by an 
MRI system. Fig. 2 is shown comprising spatial filter bank 210, decimated gradient 
phase encoding 220, interpolating system 230 and synthesis filter bank 240. Each 
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component is described in further detail below. It is to be noted that the array of RF 
detector coils forms the spatial filter bank. 

[0024] Spatial filter bank 21 0 is formed with the respective sensitivity profiles for 
spatially filtering the plurality of detected MR signals. In an embodiment, the 
sensitivity profiles are represented by F(z), where F(z) equals: 



F 2 (z) 



Equation (1) 



[0025] 
[0026] 



where F (z), ... , F (z) are the sensitivity profiles of coil 1 to coil N. 
1 N 

The overall detected signal S'(z) after the spatial filter bank F(z) can be 

Equation (2) 



described as: S'(z) = F(z)S(z) = 



[0027] where S1(z), ... , SN(z) are the original signals, and S ' l(z), ... , S ' N(z) are the 
signals after spatially filtered by sensitivity profiles of each coil. Desirably, the 
detector coils of the array are arranged to optimize the spatial encoding of the spatial 
filter bank. 

[0028] Decimated gradient phase encoding 220 collects a reduced set of phase-encoded 
MR signals resulting in fewer lines of k-space. Interpolating system 230 interpolates 
the plurality of RF and gradient-encoded signals to generate a plurality of interpolated 
signals. Alternatively, the decimated gradient encoding comprises collecting sparse 
subsets of non-rectilinear trajectories in k-space, for example fewer interleaves of an 
interleaved spiral trajectory or fewer radial lines of a radial trajectory. The signal 
generated after RF and gradient encoding and interpolation is represented by: 



S"(z) = — YSXzW m ) =— ^FizW^SizW") Equation (3) 



[0029] where M is the decimation factor (which is less than or equal to the number of 
coils or elements in the array), and W=exp(-j2 tt /M). 
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[0030] The image can be reconstructed by applying a lapped transform or filter bank to 
reconstruct interpolated signals of Equation (3). The description is continued first with 
reference to the filter bank method of reconstructing images and then with reference 
to the lapped-transform method. 

[0031] Synthesis filter bank 240 transforms the interpolated signal to generate a 
corresponding image. The signal generated by synthesis filter bank 240 is 
represented by: 



S(z) = G(zf S n (z) = -7 £G(z) T F(zW m )S(zW m ) 



=-^-(6,(0 G 2 (z) - G N (z))F 0 (z) 



S(zW) 



Equation (4) 



[0032] 



where G (z), ... , C (z) are the designed synthesis filters, superscript T denotes 
1 N 

conjugate transposition , 

su) 

represents the signals after parallel reconstruction, and the aliasing-component 
matrix is 
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F.(z) = 



(F 0 (z) F 0 (zW) F 0 (zW w - 1 )^ 
F,(z) F t {zW) - FKzW"- 1 ) 

^(z) F„(zW) - F J ,(zW M - , )J 



Equation (5) 



In an embodiment, the image S(z)is substantially free of abasing and 
amplitude distortion. To obtain S(z) free of aliasing, both spatial filter bank 210 and 
synthesis filter bank 240 should satisfy 

G(z) T F a (z) = (C(z) 0 . 0), Equation (6) 

where C(z) is an arbitrary filter function. However Equation (6) ensures only 
that there is no aliasing in the final signal, but does not guarantee absence of 
amplitude distortion in S(z). The condition for an image which is free of both 
aliasing and amplitude distortion is 

G(z) T F a (z) = {z- d 0 0), Equation (7) 

where z^ is an arbitrary delay. 

The above equation can be realized by derivation of a synthetic filter bank 240 
based on the aliasing component matrix, 

G(z) T ={z~ d 0 ■ 0)f; l (z) Equation (8) 



[0033] However if there are any poles in 

f;'(z) 

, the reconstruction system may not be stable. 

[0034] 



To guarantee a stable perfect reconstruction (as used herein "perfect 
reconstruction" refers to the condition where the signals after reconstruction have the 
same amplitudes and phase as the original signal with only a time delay), an 
intermediate filter bank V{z) (not shown) can be inserted between the decimated 
gradient phase-encoding block 220 and the interpolation block 230 for stabilizing the 
reconstruction. Thus Equation (3) becomes 
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From Equation (4), (5) and (6), the perfect reconstruction condition will be 
G(z) T V(z M )F a (z) = {z- d 0 - 0) Equation(lO) 
Therefore the synthetic filter bank 240 is represented by: 
G(z) T = {z~ d 0 ». 6]f; 1 (z)V' 1 (z m ) Equation (11) 

[0035] where V(z) can be designed to keep G(z) stable. 

[0036] In many cases of parallel MRI, especially for arrays containing a large number of 

detector elements, the sensitivity profiles of the array are overlapping. In such cases G 
(z) can be derived much more easily from the lapped transform. The manner in which 
lapped transforms may be used for deriving G(z) is described below with reference to 
Fig. 3. 

[0037] The lapped transform is a filter bank in which the impulse responses of the 

synthesis filters are the LT basis functions, and the impulse responses of the spatial 
filters are the time-reversed basis functions. Fig. 3 is an exemplary embodiment of a 
lapped transform filter bank. 

[0038] 



In FIG. 3, P is represented by 
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p, 



p = 



Po 



Equation (12) 



PJ 

Po is a LxM matrix, L=DM, and D is the overlap factor. In Pi and P2, L=(D-1)M. 
The relation between input and output signals is given by 

x(*) = PP T x(*) Equation (13) 

The superscript T denotes complex-conjugate transposition. The perfect 
reconstruction condition for the lapped transform from Equation (13) is: 
PP T =I Equation (14) 

Equation (14) implies that the basis functions of PO (the columns of PO) not only must 
be orthogonal to one another, 

P 0 P 0 T =I, Equation (15) 

but also must be orthogonal to the overlapping portion of the basis functions from 
neighboring blocks, 

P 0 aP o T =0, Equation (16) 

where the one-block shift matrix is 



(O = ^ * j Equation (17) 

and the identity matrix I is of order (D-l)M. 



[0039] In a further embodiment, the detector elements of the array (RF coil assembly 1 47 
of Figure 1) comprise linear strips. Such an embodiment is referred to hereinafter as 
SAFARI or Strip Array For Analytic Reconstruction of Images. As used herein, a strip 
array refers to a coil array in which the conductive elements are microstrips of 
conductive (e.g. copper) material of a selected length connected to a ground plane 
and arranged to form a planar or cylindrical shape. In general, there are 8 or 16 strips, 
but it is to be appreciated that variations known to one skilled in the art could also be 
employed. 

[0040] 



Inherent with strip arrays , the signal from each strip introduces a different 
frequency modulation, as shown in Fig. 4. The mathematical representation of Fig.4 is 
illustrated in Fig. 5. 
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From Fig. 5, it is seen that: 

p, (m) = f{mM + 0 Equation (18) 

q i (m) = g(mM + i) * 

[0041] Where p(m), and q(m) represent mathematical equivalents for spatial filter bank f 
(mM + l) and synthesis filter bank g(mM+l), respectively. Here i+0, 1 , M-l . 

[0042] When a large number of strips are concentrated in a region that is relatively small 
compared to the sensitivity profiles, Fig. 5 can be simplified to Fig. 6, where the 
parallel reconstruction is simply a block-in From Fig. 5, it is seen that:verse discrete 
Fourier transform (DFT). 

[0043] When the RF coil assembly is a strip array, each strip has an explicit phase 

relationship relative to a spatial location within the object being imaged. Referring to 

Fig. 4, the relationship is expressed as W kj =exp(j2 tt k/M) for M=0 to (M-l), where 

M 

M is the number of strips in the array. Referring to Figures 4-6, the explicit phase 
relationship in the sensitivity profiles can be used for simplifying the lapped transform 
reconstruction described earlier. 

[0044] In further embodiments of the present invention, reconstruction may be improved 
by adjusting the sensitivity profiles of the coils of the RF detector array (RF coil 
assembly 147 of Figure 1) to realize an orthogonal sensitivity profile. Orthogonality is 
useful for reducing the matrix condition number, thereby making the reconstruction 
more robust. In a first embodiment, the RF detector array is a planar strip array 
comprised of multiple conductive strips and coupling capacitors interconnecting the 
strips. By adjustment of the interconnecting capacitors and input ports for each 
respective strip, the desired sensitivity profile can be obtained. In a second 
embodiment for adjusting sensitivity profiles, post processing can be used to improve 
orthogonality. 

[0045] While only certain features of the invention have been illustrated and described 
herein, many modifications and changes will occur to those skilled in the art. It is 
therefore to be understood that the appended claims are intended to cover all such 
modifications and changes as fall within the true spirit of the invention. 
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